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MONOTONE SOBOLEV MAPPINGS 

of planar domains and surfaces 

TADEUSZ IWANIEC AND JANI ONNINEN 


Abstract. An approximation theorem of Youngs (1948) asserts that a con¬ 
tinuous map between compact oriented topological 2-manifolds (surfaces) is 
monotone if and only if it is a uniform limit of homeomorphisms. Analogous 
approximation of Sobolev mappings is at the very heart of Geometric Function 
Theory (GET) and Nonlinear Elasticity (NE). In both theories the mappings 
in question arise naturally as weak limits of energy-minimizing sequences of 
homeomorphisms. As a result of this, the energy-minimal mappings turn out to 
be monotone. In the present paper we show that, conversely, monotone map¬ 
pings in the Sobolev space , 1 < p < oo , are none other than -weak 

(also strong) limits of homeomorphisms. In fact, these are limits of diffeomor- 
phisms. By way of illustration, we establish the existence of energy-minimal 
deformations within the class of Sobolev monotone mappings for p-harmonic 
type energy integrals. 


1. Introduction 

There has been recently increasing interest in the Sobolev - homeomor¬ 

phisms and their weak and strong limits. In planar domains (or surfaces) and 
p ^ 2, these limits turn out to be monotone mappings in the topological sense 
of C.B. Morrey [42], see Definition l.I below. We shall see that, conversely, every 
^i:P-map that is continuous and monotone between Lipschitz domain (this time 
for any exponent 1 < p < oo ) can be approximated by homeomorphisms uniformly 
and strongly in the Sobolev norm; hence, by '^“-diffeomorphisms, also by piecewise 
affine homeomorphisms. A motivation for Sobolev mappings comes from the study 
of extremal problems in Geometric Function Theory (GFT) [4, 20, 22, 25, 27, 29, 30]. 
Further motivation comes from the fields of materials science such as Nonlinear 
Elasticity (NE) [2, 5, II, 39, 47, 48]. 

Let us begin with a brief analysis of the p-harmonic energy of homeomorphisms 
/ : X Y between bounded domains X, Y C . 

(l.I) <ffp[f]= f \Df{x)\Pdx , where \Df{x)\^ ^Tr[D*fix) D fix)] 

Hereafter Dfix) stands for the Jacobian matrix of / (deformation gradient). De¬ 
note by c^(X, Y) the class of orientation preserving homeomorphsms / : X Y 
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of finite energy and 

( 1 , 2 ) 

The infimum may or may not be attained. If not, we wish to find a map ho G 
Y) in close proximity to J^(X, Y) such that 

(1.3) J \Dho{x)\P dx = S’p(X,Y) 

It is natural to look for ho as -weak limit of an energy-minimizing sequence 
of homeomorphisms hj : X Y. That this indeed would solve the problem (at 
least for Lipschitz domains and p ^ 2) is far from being obvious. We always have 
the upper bound 

(1.4) J \Dho{x)\P dx <^p(X,Y) 

Equality occurs if and only if the minimizing sequence hj ho actually con¬ 
verges strongly in , which places ho in the closest possible proximity 

to J^p{X,Y). Equality (1.3) is of practical significance. Indeed, once ho fails to be 
injective, it will tell us when to stop the minimizing sequence of homeomorphisms; 
that is, prior to the conditions favorable for the collapse of injectivity. This is a 
phenomenon known as interpenetration of matter. 

Monotone Mappings. The concept is due to C.B. Morrey [42] (1935). The 
interested reader is referred to the Proceedings of the Conference on Monotone and 
Open Mappings [40] (1970) and the series of early papers by G.T. Whyburn for 
further reading. 

Definition 1.1. Let A and B be compact metric spaces. A continuous map 
h: A B is said to be monotone if every fiber h~^{b} C A of a point 6 G B is 
connected. As shown by G.T. Whyburn the preimage h~^{C) = {a G A: h{a) G C} 
of any connected set C C B is connected in A . 

Let us emphasize that monotone mappings are continuous, by the definition. 

While manifolds and Riemannian metric tensors are not the primary issues, it 
is desirable to keep them in mind since the topological aspects really crystalize in 
the manifold setting. Thus we choose and fix, as reference manifolds, two - 
smooth closed (compact without boundary) oriented Riemannian 2-manifolds 3^ 
and 3 of the same topological type. We shall consider multiply connected Jordan 
domains X C and Y C 3. Precisely, X and Y will be obtained by removing 
from 3" and 3 the same number, say 0 < £ < oo, of closed disjoint topological 
disks. In fact, any pair of topologically equivalent open -smooth surfaces with 
i boundaries can be obtained in this way. Such are planar multiply connected 
Jordan domains in ~ C . 

Theorem 1.2 (Youngs, [51]). Let X C 3 and Y C 3 be Jordan domains of 
the same topological type. A map h : X Y is monotone if and only if it is a 

uniform limit of homeomorphisms hj : X Y ,j = 1, 2,.... 

It follows, in particular, that the boundary mapping h : dX dY is monotone 
as well. 
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Sobolev Variant of Youngs’ Approximation Theorem. The main result: 

Theorem 1.3. Let H C and Y C be Jordan domains of the same topological 
type, Y being Lipschitz. For every monotone map h : X Y in the Sobolev 
space W^'P{1L,Y),1 < p < oo, there exists a sequence of monotone mappings 
hj : X Y such that: 

(i) hj : X Y are homeomorphisms 

(ii) hj h uniformly on X 

(iii) hj —>■ h strongly in Y) 

(iv) hj = h dY, for j = 1 , 2, ... 

Let (X, Y) denote the class of orientation preserving monotone mappings 

/ : X Y of finite p-harmonic energy, 1 < p < oo . Theorem 1.3 implies (by 

direct method) that there always exists h G ^p (X, Y) with smallest p-harmonic 
energy. More importantly, the energy of h equals precisely the infimum of the 
energy among homeomorphisms: 


f \Dh{x)\Pdx= min_ / \Df{x)\Pdx= inf / |L>/(x)|^’dx = 4 (X, Y) 
Jx f€^p{X,Y)Jx Jx 

see Section 6 for a definition of the p-harmonic integrals on surfaces. In other words, 
no Lavrentiev Phenomenon occurs in the class of monotone Sobolev mappings. 

Remark 1.4. It is worth noting that, for an arbitrary pair (X, Y) of topologically 
equivalent planar domains, diffeomorphisms are dense in Jtfp(K,Y ), see [24]. Thus 
one can take for hj in Theorem 1.3 a sequence of “^“-smooth diffeomorphisms. 
Furthermore in case p = 2 , if one is willing to sacrifice the boundary condition (iv) 
then the diffeomorphisms hj : X Y can be chosen to be homeomorphisms up 
to the boundary, again denoted by hj : X Y, see [23]. 

Remark 1.5. It is not difficult to see that, in case of planar Lipshitz domains, a 
sequence of homeomorphisms hj : X Y which converges weakly in Y) , 

p ^ 2 , actually converges uniformly to a monotone mapping h : X Y. Thus, 
in particular. Theorem 1.3 implies that h can be realized as #'^’^’-strong limit of 
homeomorphisms. We therefore recover a result in [31]; accordingly, weak sequential 
closure and strong closure of J^p{X,Y ), p ^ 2, are the same. 

2. Topological Preliminaries 

This section is intended as a gentle introduction to underlying geometric analysis 
on planar domains and surfaces. When discussing the Sobolev class W^'P{X,Y) , 
it is particularly convenient to embed the reference manifolds and '3^ into 
an Euclidean space. We may, for example, use the Nash-Kuiper embedding the¬ 
orem [36, 43] which assures that the oriented 2-manifolds are isometrically 
-embeddable in . Thus we may (and do) assume that 

SJ ,3^ C ; inclusion being a - embedding. 

Note that for the purpose of defining #'^’^’(X, Y) one needs only assume uniform 
upper and lower bounds on the metric tensors. 
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When interpreting the conclusions, one might view the surfaces X and Y as thin 
films in , or flat plates in case = ^3^ = ~ §2 ^ 

The components of ^ \ X and \ Y are topological disks, say 

^ \ X = Xi U X 2 U ... U X<> and ^ \ Y = Yi U Y 2 U ... U Y^ 

d0f d0f 

Their boundaries, denoted by &%.i, = and 9Yy = , are exactly the com¬ 

ponents of 9X and dY , 

ax = Xi U ... U Xi and SY = Ti U ... U Te . 

By convention, X and Y have no boundary if ^ = Q, uy which case 'K = 3^ and 
Y = 

Continuous Functions and Homeomorphisms. We shall work with various 
function spaces defined on subsets of the reference manifold 33 : 

• For a compact subset Kd 33 vie denote by ^ (A) the space of continuous 
functions h: A —>■ R^ furnished with the norm: 

II *lk(A) = max|h(a;)| 

x£A 

The notation h G "^(A, B) will be used if we want to make explicit the 
range of h : A —>• B C R^ . 

• J^(X, Y) consists of orientation preserving homeomorphism h: X Y. 
Every homeomorphism h : X Y gives rise to a one-to-one correspon¬ 
dence between boundary components Xi, X 2 ,..., Xe and Ti, T 2 ,..., Te by 
means of cluster limits. We conveniently rearrange the indices so that the 
boundary correspondence reads as follows: 

(2.1) h:X^-^T,y , for v = 1,i ] 

This arrangement will tacitly be assumed throughout this paper for home¬ 
omorphisms in the class J^(X, Y) and their uniform limits. 

• J^(X, Y) consists of homeomorphisms h : X Y which extend contin¬ 

uously to the closure of X. Note that the continuous extensions become 
monotone maps from X onto Y, still denoted by h: X Y. They 
take 5X onto dY ; specifically, h{X3) = for v = 1, ...,£. We refer to 
h: X^ , II = 1, ...,£, as the boundary maps. Each of these bound¬ 

ary maps is monotone. 

• Jf’(X,Y) is the space of homeomorphisms h: X Y. 

Similar notation, with the obvious analogous meaning, will be used for spaces of 
mappings defined on other subsets of 33 . 

Monotone Mappings Versus Uniform Limits of Homeomorphisms. Let 

Kd 3333 and B c be compact. The class 

./#(A, B) C (A, B) stands for the space of monotone mappings h: A B . 
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Obviously M’iKM) C B) • We now invoke the Approximation Theorem 

of J. W. T. Youngs [51]. 

Theorem 2.1. Let X and Y he topologically equivalent compact 2-manifolds. 

Then for every monotone map h: X Y there exists a sequence of homeo- 
morphisms hj : X Y converging uniformly to h. In symbols, 

..#(X,Y) = 

This theorem will provide us with a powerful tool when dealing with monotone 
mappings. We shall appel to it repeatedly to either recover or refine the well known 
properties of monotone mappings betwen Jordan domains 'K G 2T and Y G . 
For example, Theorem 2.1 readily implies that 

Lemma 2.2. We have the inclusions 

Jf{%Y) 5 ^^(X,^) £ .^(X,Y) = .^(X,Y) g ^{%Y) 

Moreover, for h G ./#(X, Y), each boundary map h : dXi, i9T^ is monotone. 

We shall also take advantage of the following refinement of the Modification 
Theorem in [51]. 

Lemma 2.3 (Homeomorphic Extension of a Monotone Boundary Map). Let Xo C 
and Yo G'SL be simply connected Jordan domains. Then every monotone map 
h : 9Xo dYo admits a continuous monotone extension h : Xo^^^ Yq. Such 
an extension can be further modified to become a homeomorphism between Xo and 
Yo. 

Proof. First, with the aid of the uniformization theorem, we transform Xo and Yo 
homeomorphically onto the closed Euclidean disks. We obtain a monotone map 
between circles, which we extend (in a radial fashion) to a monotone map of the 
disks. We then lift such an extension back to the reference manifolds ^ and , 
completing the proof of the first statement. Then, by ’’Modification Theorem” in 
[51], this extension can be modified to become a homeomorphism between Xo and 
Yo. □ 

Remark 2.4. In Theorem 2.1, if one is willing to forgo the univalence of the boundary 
mappings hj : 9X dY , but only wants them to be injective from X Y, 
then we can ensure that hj = h on 9X. However, this variant of Youngs’ theorem 
will not be exercised here. 

Monotone Extension Outside X . 

Lemma 2.5. Every h G .^(X, Y) extends as a continuous monotone map between 
reference manifolds. 

Proof. If X has no boundary then so does Y. Hence X = and Y = ‘3^. 
Now, consider a boundary component Xi, G 5X and the corresponding boundary 
component Tj, C dY. Recall that h: Xi, is monotone. Both X^ and 

are topological circles on the reference surfaces and 31, respectively. Filling 
these circles with topological open disks, say Xj, C and Y,, G 31, results in 
closed 2-cells in the reference manifolds. The extension is immediate from Lemma 
2.3. □ 
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Remark 2.6. With little extra efforts one can ensure that such an extension is a 
homeomorphism between the components X,y and Yi ,, for all v = 1,2, . But 

there will be no need for this. 

We continue to use the same notation h : ^ ‘3^ for the monotone extension. 

Thus h becomes an element of ,3^) as well. There will be no need of any 

regularity ol h ■. 3^ outside X . 

Cells and Pre-cells. A cell in the reference manifold 3^ is any simply connected 
Jordan domain Q <Z 3^. Its preimage h~^{Q) C 311 under a monotone map 
h : 31 3 is still simply connected but not necessarily a Jordan domain. We 

refer to h~^{Q) CL 31 && pre-cell in 31. Beware that not every simply connected 
domain in 31 comes as a pre-cell for h ; and that, in general, the closure of a 
pre-cell may not contain h~^{Q). In fact, the strict inclusions h~^{Q) h~^{Q) 

and dh~^{Q) ^ h~^{dQ) are typical of monotone mappings. Associated with 
C 31, Y C 3 and the monotone map h : X Y are the concepts of internal 
and boundary cells. 

Definition 2.7. The term cell in a domain Y C 3 refers to any simply connected 
Jordan domain Q C Y which satisfies one of the following conditions: 

• Q (£ Y, we call it an internal cell. 

• Q n dY = C is a closed Jordan arc (not a single point). We refer to such 
Q as boundary cell in Y and to C as its external face. In this case we shall 
also make use of the set Q-|. = Q U C . 

Every boundary cell Q C Y can be extended beyond its external face to be¬ 
come a cell in the reference manifold 3. Let Q* C be such an extension so 
Q= Q* n Y. Here dY splits Q* into two simply connected Jordan domains. 

To every cell Q in Y (internal or boundary) there corresponds a simply con¬ 
nected domain in X, called pre-cell in X. It is defined by the following rule. 

• If Q (£ Y, then we call 14 = h~^{Q) (£ X the internal pre-cell in X. 

• If Q n i9Y ^ 0 , then we call U = h~^{Q*) n X the boundary pre-cell in 
X . It is of no importance which extension Q* is taken in this formula for 
14. In fact, we always have 14 = h~^{Q+) fl X. 

The reader is cautioned that in general lA is not the preimage of the boundary cell 
Q C Y; we have only the inclusion ft.“^(Q) C U. Furthermore, the pre-cells in X 
need not be Jordan domains, and that is why some complications (topological and 
analytical) are to be expected. 

Lemma 2.8. The pre-cells (both internal and boundary) are simply connected do¬ 
mains in X. 

Proof. We again take advantage of Youngs’ Approximation Theorem. The lemma 
holds a h : 31 is a homeomorphism. But it also holds if h is monotone; 

just use the sequence of homeomorphisms converging uniformly to ft., to represent 
the pre-cell as union of an increasing sequence of simply connected domains. □ 

Our next topological fact, which actually strengthens Lemma 2.3, will require 
some work. Recall that we have extended ft G ./#(X, Y) to a map in 3^{31,3). 
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Lemma 2.9 (Homeomorphic Replacement in Pre-cells). Let Q C Y be a cell in 
Y (internal or boundary) and W C X the corresponding pre-cell in X . Then there 
exists a monotone map G such that 

• h^ :U Q is a homeomorphism 

• h^ = h ■. SL\U 

Proof. The proof comes down to monotone mappings between 2-spheres and a re¬ 
sult by T. Rado [46] page 66, 11.1.47. Let us consider two cases. 

Case 1 [internal pre-cell]. Suppose lA = h~^{Q) where Q g Y. Choose and 

fix slightly larger cell Q (s Q' (s Y, which gives us a larger pre-cell U U' = 
h~^{Q') (£ X. We view U' and Q' as open simply connected Riemann surfaces. 
Note that the map h : U' ^ Q' is continuous and takes dW into dQ'. 

Denote by 14' and Q' the Alexandroff one-point compactifications of W and 
Q', respectively . These are topological 2-spheres. The unique continuous exten¬ 
sion h : W Q' of h : U' Q' remains monotone. At this point we 

may appeal to [46] page 66, II.1.47, which asserts that there is a monotone map 
h^^ : 14' Q' which takes U homeomorphically onto Q and agrees with h 
outside U. 

Case 2 [boundary pre-cell]. Suppose U = h~^{Q*) DU, where Q* is a cell 

H p'F 

in '3^ such that Q = Q* fl Y and Q* fl dY = C is an open Jordan arc. Denote 
by D = h~^{Q*) the pre-cell in under the map h : SL 3^. As in the 

previous case we find a monotone mapping H = ho, : 4^ 3^ that agrees with 

h : lAP \Ll 3^ \ Q* and takes D* homeomorphically onto Q* . The issue is 

that H : D Q* need not take W = D fl X onto Q* fl Y and, even if it does, 

need not coincide with h outside U. The idea is to correct H within D. For, 
we look closely at the crosscut of D by the boundary of X, say by the component 
X = Xy C 9X, for some v = !,...,£. Let T = denote the corresponding 
boundary component of dY. Recall that the boundary map h : X T is also 
monotone; that is, the preimages of connected sets in T are connected in X. In 
this way we have defined an open subarc of X, 

r = {cc G X; h{x) G C }. Note that T X, since C = h{X). 

This subarc has two endpoints (the limit points), say a, & G dLl. It is a topological 
folklore that an open Jordan arc in a simply connected domain D, whose endpoints 
lie dLl, splits D into two simply connected subdomains. These subdomains are: 

W = D n X and V = D \ X , hence a decomposition D = W U T U V 

Homeomorphism H : D Q* yields a decomposition of the cell Q* C 3^; 
namely, 

Q* =U^UT^UV^ , where = H{U) , T^ = H(T) and V„ = H{V) 

The open Jordan arcs C <Z Q* and T^ C Q* share common endpoints which we 
denote by A = h{a) = H(a) G dQ* and B = h{b) = H{b) G dQ* . We note that 
Uu and Vu are Jordan domains for which F^ U {A, B} constitutes their common 
boundary. Precisely, we have 

dU^=aUT^ , av^=/3UF^ , dU^ndV^ = T^ 
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where a is the closed sub-arc of dQ* between A and B that lies in Y, and /3 is 
the closed sub-arc of dQ* between A and B that lies in \ Y. 

Now, just the fact that h : H Q* is continuous and agrees with H on dO, lets us 
observe that the mapping ho : Q* Q* (not necessarily injective) extends 
continuously as a monotone map of Q* onto itself. Upon such an extension, the 
boundary map h o H~^ : dQ* dQ* becomes the identity. Let us see how this 
extension acts on dlA„ . It is still the identity map on a C dQ* and it takes the 
closed subarc L^C dUjj monotonically onto C C dQ* . Thus we have a monotone 
map ho H~^ : dU^j d{Q* fl Y). It is important to observe that both U„ and 
Q* n Y are Jordan domains. At this stage we appeal to Lemma 2.3 . Accordingly, 
we extend h o H~^ : dUg d{Q* fl Y) continuously, and as a homeomorphism 
inside the curves. Denote the extension by [h o H~^)^ : Ug Q* fl Y. In 
summary, we have constructed a continuous monotone map F : ^ ^, 

( identity in ^ \ Q* 

U = < hoH~^ AiFg 

{ (hoiL-i)tt in'^^ 

The composition hg = F o H : ^ is the desired replacement, as claimed 

in Lemma 2.9. □ 


3. Analytical Requisites 

Since the reference manifolds 3^ and are closed oriented Riemannian 2- 
manifolds of class , we may speak of the Sobolev class of mappings 

h : V, ^ iF defined on any open subset C 3^ . We do not reserve any particular 
notation of the metric tensors on 3" and '3 , though we fix them for the rest of 
this paper. The volume element on 3 , denoted by dx , is the one induced by the 
metric tensor. We recall the C^-isometric embeddings 3 and 3 . 

3.1. Sobolev Mappings Between Surfaces. The Sobolev space of 

real-valued functions on D, (Z 3" will be endowed with the seminorm 

\\<l>\\wnp{Q,) = (^J^\D(l){x)\P dx'^ , l^p<oo 

where |I2^(a;)| , defined almost everywhere, stands for the norm of the linear tangent 
map D(j){x) : T2,(r2) —>■ K. with respect to the inner product in T 2 ,(D). The Sobolev 
class C consists of mappings h = (h^,: D —>• 

whose coordinate functions belong to W^''P{Vt) and h{x) G 3 for almost 

every a; G S2. It may be worth reminding the reader that for 1 ^ p < 2 = dim Q 
topology can inhibit the space from being dense in W^’'‘’{il.,3). The 

interested reader is referred to [10, 17, 18, 16] for this issue. This problem, however, 
disappears completely since our mappings in question are continuous. De facto, 
when 1 < p < 2 , the continuity assumption of monotone Sobolev mappings is 
critical for the subsequent arguments; it is superfluous, however, when p ^ 2. 
Royden p-Algebra. The class 3P{fl) = ‘^(D) fl consist of real-valued 

functions in the Sobolev space that are continuous on O. This is a 

Banach algebra with respect to the sub-multiplicative norm 

lll^ll|5?p(n) = ll'/'ll<g’(n) + ll^^ll5rLp(n), 11^^'• V'llli*p(n) ^ lll^^llli*p(n) ■ IIIV’lll^*p(n) 
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3.2. p-Harmonic Boundary-Value Problem. We record less familiar aspects 
of the Dirichlet problem for the p-harmonic equation. 

(3.1) div= 0 , for (/) G , 1 < p < oo 

in planar simply connected domains. In general, simply connected domains may 
have rather odd boundary. 

Lemma 3.1. Let O C be bounded simply connected domain and $ G . 

Then there exists unique </> G that is p-harmonic in LI, 1 < p < oo, and 

agrees with $ on dil. If ^ G , then also f G . Moreover, 

(/) G $ + and [ \V(t)\P < [ |V4>|P 

JQ Jfl 

Equality occurs if and only if ^ f. 

Proof. We will only briefly outline the key points of the proof. For a thorough 
treatment of the Dirichlet problem we refer the reader to [19] and [38]. By virtue 
of the Wiener’s criterion the first statement of the lemma holds whenever the 
complement \ D is p-thick at every boundary point, see [38, Corollary 6.22] 
and [38, (2.22)] for a formulation of Wiener’s criterion. Simply connected domains 
indeed satisfy this criterion, a fact not difficult to verify, but it is not explicitly 
exemplified in the vast literature. For the second statement we minimize the p- 
harmonic energy in the class $ + (so-called variational formulation) which 

does not require any regularity assumption on the domain Ll. The only point is 
to show that the variational solution extends continuously to the boundary and 
it coincides with $. This again follows by Wiener’s criterion, see [19, Theorem 
6.27]. □ 

3.3. p-Harmonic Replacements. Let a map F = t6-|-fv:f2—defined 
in a domain D C C, belong to the Sobolev space (O, C) , 1 < p < oo . 

Definition 3.2. F is said to be p-harmonic (coordinate-wise) if both coordinate 
functions u and v satisfy the equation (3.1). We shall introduce the unisotropic 
p -harmonic energy of F: 

(3.2) ^[F] = [F] ^ (lVuJJ’ + iVnl^’) 

and repeatedly abbreviate this as energy of F ; because the exponent p will remain 
fixed. 

Remark 3.3. This terminology is different from what can be found in the literature; 
the term p-harmonic mapping is usually reserved for the coupled p-harmonic sys¬ 
tem divjZl/j^’“^Il/ = 0. There is a subtle distinction between these two concepts. 

Recall from Lemma 3.1 that for any F G .^p( 0,]R^) in a bounded simply con¬ 
nected domain its boundary map F : dLl —>■ admits unique p-harmonic exten¬ 

sion to D . The question arises whether such an extension is injective. The answer 
depends on how F runs along dVl. Of course, the boundary map F : dTl —>■ R^ 
must admit at least one homeomorphic extension inside D. This is also sufficient 
if F(dLt) is a convex curve. 
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Lemma 3.4. Let C be bounded simply connected domain and let F G 
be p- harmonic (cordinate wise) in ft. Suppose there is 4* G '^($7,®.^) 
that agrees with F on dfl and takes homeomorphically onto a convex domain 
A . Then F is a -diffeomorphism of O onto A . Moreover, 

(3.3) £’n[F] < , whenever <S> G F + C). 

In this setting 4* plays the role of the classical Dirichlet boundary data. It helps 
to capitalize on the topological properties of F near the (rather weird) boundary 
dLl. We emphasize that 4/ is not required to have any Sobolev regularity in LI. 


Remark 3.5. Lemma 3.4 is a p-harmonic analogue of the celebrated Rado-Kneser- 
Choquet Theorem [14]. Under additional assumptions on the domain 11, a p 
-harmonic (coordinate wise) analogue of Rado-Kneser-Choquet Theorem was first 
shown by Alessandrini and Sigalotti [1], see [21] for the isotropic case. The paper [1] 
is concerned with the domains which satisfy the external con condition. This would 
be redundant, as mentioned in ’’Remark 3.2” of this paper. However, the essential 
shortcoming is that neither formal statement nor the proof in case of non-Jordan 
domains and non homeomorphic boundary data are provided in [I]. Thus we must 
work out additional arguments. 


Proof. We may assume, in addition to the hypotheses above, that 41 takes LI 
diffeomorphically onto A . This is permissible by a theorem of Rado [45] , see 
also [41], which asserts that to every homeomorphism 4i : H — >■ and continuous 
function r : H —>■ (0, oo) there corresponds a diffeomorphism 417 . : H —>■ such 
that |4'(x) — 4 't-(x) I < t{x) in LI. Therefore, one can replace 41 in Lemma 3.4 by 
a diffeomorphism 4ir with t{x) = dist(a;,9H). 

Now, consider an increasing sequence of smooth convex domains Ai (s A 2 ^ ... ^ 
A„ (s A„+i.... (s A whose union is A. The corresponding preimages under 4i 
are smooth Jordan domains in Ll , 

H„ = 4'-i(A„), Hi (E H 2 g ... g (£ H„+i...., |Jh„=H. 


Next, we replace each diffeomorphism 4i : A„ by a homeomorphism 4i„ : 

A„ which is p-harmonic (coordinate-wise) in H„ (thus a diffeomorphism in 
Lin ) and agrees with 4i on , see Theorem 5.1 in [1]. We just constructed a 
sequence of continuous mappings Fn : Ll A which are homeomorphisms on Ll 
and coincide with F on dLl , 


Fn 


41 in H \ Lin 

41 n in Lin 


, with the Jacobian determinants J{x, 4i„) > 0 in 


This sequence converges uniformly to F . Indeed, for every x G Ll we have 
(3.4) \Fn{x) - F{x)\ ^ V2 ||4i - ^ 0 

The latter estimate is trivial when a; G H \ H„ . To see that this also holds for 
X G Lin we argue by a comparison principle [19, Comparison principle 7.6, page 
133] in a straightforward way. Namely, the coordinate functions oi F = {u, v) and 
Fn = {un,Vn) , being p-harmonic in H„ , satisfy 

\u„{x) - u{x)\ < \\un - u|LfnQ a^d \Vn{x) - v(x)\ < ||t>„ - ll|L/ao -i 
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which yields the desired estimate in (3.4). Concerning the positive sign of the 
Jacobian determinant of F, we shall appeal to a p-harmonic variant of Hurwitz 
Theorem [31, Theorem 4.9]; that is, 

Theorem 3.6. If a sequence —>■ K.^ of p-harmonic (coordinate-wise) orien¬ 

tation preserving diffeomorphisms converges uniformly to F : H ^ "Mf in a domain 
n C , then either J{x,F) > 0 everywhere in ft or J(x, F) = 0 in . 


Let us first exclude a possibility that J{x,F) = 0 in . For this, we choose and 
fix a nonnegative test function rj G (A) whose integral mean equals 1. Let 
G (s A denote the support of p. Since Fn ^ F (uniformly) and F{dlt) = dA it 
follows that 


^(G) (s G . In particular, ^(G) g g G , for sufficiently large fc 

n^l n^l 

Hence, for n > k, we have 



p{Fn{x))J{x,Fn)dx 


iFr^ink) 


v{y)dy ^ 


p{y) dy =1 . 


Recall that Fn ^ F on G^+i gj . Since Fn are p-harmonic in G^+i , we also 
have DFn =4 DF on . This follows from local “^^’“-estimates of p-harmonic 
functions [26, 49, 50]. It is now legitimate to pass to the limit in the above estimate 
to obtain J^^p{F{x))J{x, F) dx > 1. In particular, J{x,F) ^0 in G^ . Thus 
J{x, F) > 0 everywhere in G^ . But k can be as large as we wish, so J{x, F) > 0 
everywhere in G. 

In particular, F is a local diffeomorphism in G. On the other hand the map 
F : G A is a uniform limit of homeomorphisms. Therefore, F is a global 
diffeomorphism, completing the proof of Lemma 3.4. □ 


4. Approximation inside a given pre-cell 

From now on Y is a Lipschitz domain. This means that every point po G dY 
has a neighborhood O C ^ and a local (“^^ -smooth) chart k : O —>■ such 

that k{0 n dY) is a graph of a Lipschitz function. Regarding the given map 
h G “^(X, Y) n 1^^’P(X, R^) , I < p < oo, we recall that it extends as a monotone 
map h : SF ^. No regularity outside X is required, as this extension will 
assist us only in the topological aspects of the proof. 


As a preliminary step in the proof of Theorem 1.3 we shall approximate h : 
X Y with Sobolev mappings which are univalent within a given pre-cell, and 
remain unchanged outside the pre-cell. The approximation is understood by means 
of the metric in the Royden space I%p(X,Y) = '^(X,Y) fl 1^^’P(X,R^); that is, 
uniformly and strongly in 1^^’^’(X, R^). Precisely, we have: 

Proposition 4.1. Let Q C Y be a cell in Y, which we assume to be Lipschitz 
regular, and lA be its pre-cell in X. Then there exists a sequence of monotone 
mappings hj : X Y such that: 

(a) /ij G‘r(X,^niri’P(X,R3), j = I,2,... , 

(b) hj = h on X\U , 

(c) hj : U Q are homeomorphisms, 

(d) hj ^ h uniformly in X, 
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(e) hj —>■ h strongly in R^), 

Remark 4.2. Upon the extension hj = h-.^\'K '3^\Y, we obtain monotone 

mappings between reference manifolds, again denoted by hj : ‘W . 

Proof. It will takes 5 steps to complete the proof of Proposition 4.1. 

Step 1. Reduction to simply connected Jordan domains in . We view as 
the extended complex plane C = CU {oo} ~ equipped with the chordal metric. 
Choose and fix a simply connected neighborhood Q' §) Q such that Q' fl dY 
becomes a Lipschitz regular Jordan arc. Obviously such Q'does exist. Thus Q'flY 
is a boundary cell in Y, even if Q was an internal cell. Denote by W = h~^{Q') C 
JP . Since h is monotone, U' is a simply connected neighborhood of the continuum 
h~^{Q) (E W and h : dW —>■ dQ'. Take a look at the commutative diagram 

W - h -^ Q' 

I I 

(j) if 

i ^ i 

R^ - h -R^ 

where : W R^ and ip : Q' R^ are -diffeomorphisms. For the 
existence of such diffeomorphisms one may appeal to the uniformization theorem 
[44] and [33, 34, 35]. It asserts that every simply connected Riemann surface is 
conformally equivalent to either open unit disk, the complex plane or the standard 
sphere; thus in our case, -diffeomorphic to R^ . Note that h = ip o h o (p~^ : 
R2 R2 extends continuously to the one-point compactification of the plane; 
that is, to a monotone mapping of the Riemann sphere onto itself, still denoted by 
h : C C , /i(oo) = oo . With the aid of stereographic projections we move the 
sets X n ‘p(U') and Y fl ip{Q') away from oo . Now the proof of Proposition 4.1 
reduces to the case in which 

jr = c, r = c 

X, Y C C are bounded simply connected Jordan domains, 

Y and the cell Q C Y are Lipschitz domains. 

Step 2. Further reduction. We shall need Q C Y to be convex; for instance, 
the unit square. For this, since Q C R^ is Lipschitz domain, one might try to use 
a bi-Lipschitz transformation F : R^ ^ that F : Q (0,1) x (0,1). 

Remark 4.3. Although the existence of such F poses no problem, a caution must be 
exercised. Every bi-Lipschitz map F : A B between bounded planar domains 
and its inverse F~^ : B A induce bounded (nonlinear) composition operators: 

Fj, : , by the rule F^{g) =Fog , 

Fj,-i : , by the rule Fff^{f) = F-^ o f , 

whatever the domain D C R^ is. But in general the continuity of these operators 
is questionable [15]. Fortunately, there is a satisfactory solution to this puzzle. For 
Lipschitz domains such as A = Q and B = (0,1) x (0,1) one can construct special 
bi-Lipschitz transformation F : R^ ^ F : A B, for which the induced 
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composition operators and ^ are indeed continuous. Actual construction of 
such F is presented in [31]. 

Thus we may assume that Q = (0,1) x (0,1) C Y. Furthermore, in case of a 
boundary cell, we assume that its external face Qn(9Y equals {(a;, 0), 0 ^ a; < 1} . 

Step 3. Covering by small cells. Given £ > 0, we cover Q by the family of 
overlapping open squares Qi, Q 2 , ■■■, Qn of diameter less then e in which every 
point in Q belongs to at most three of the closed squares in this family. In math¬ 
ematical terms, 

• Q= Q 1 UQ 2 U...UQJV, N = N{e) 

• diamQi <£, i = l,2,...,A^ 

• ^3, forzGQ 

• Each Qi is either internal or a boundary cell for Y. 

Construction of such a cover poses no difficulty. We now consider the corresponding 
pre-cells in Ui = h~^{Qi) C C X. 

• U = U 1 UU 2 U ...UUn , N = Nis) 

• 1 ^ (a;) < 3 , for X G W 

• Each Ui is either internal or a boundary pre-cell in X. 

It may be worth pointing out that the pre-cells Ui can remain very large in diam¬ 
eter as £ approaches zero. Typically this occurs when a continuum collapses into 
a point. 

Step 4 ■ A chain of p-harmonic replacements. We are going to construct by 
induction a chain /i /2 /at /n+i of mappings ft : C C. We 

set fi = h. For the induction step suppose we are given a monotone map ft : 
X Y of Sobolev class and its monotone extension fi : C C. 

Consider the cell Qi C Q C Y (internal or boundary) and the corresponding pre¬ 
cell Ui C U C X for the mapping fi. By Lemma 2.9, there exists a monotone 
map 'I'i : C C such that 5'i : Ui Qi is a homeomorphism and 'I'i = /i : 
<C\Ui C \ Qi. Then, with the aid of Lemmas 3.1 and 3.4, we can modify fi 
within the pre-cell Ui to obtain, 

f. ^ f fi in C\U^ 

\ p-harmonic of class fi F Wo'^{lAi) in Ui 

We emphasize (in view of Lemmas 2.9 and 3.4) that /i+i takes Ui homeomor- 
phically onto Qi, whereas under the map fi some points in Ui may collapse into 
^Q^. 

Here are the essential properties of these mappings. It may be worth reminding the 
reader that we are dealing with rather weird domains. 

• Each fi belongs to the Sobolev class #'^’^’(X, Y). This is due to Lemma 
3.1 which yields /*+! - /* G Wo'^{U^,<C). 

• The energies are nonincreasing; ^x[/i] ^ ^x[/i-i-i] ^ ^ <^x[/i] ■ 



14 


T. IWANIEC AND J. ONNINEN 


• Each /i+i is locally injective in UW 2 U... UWi (no branch points). This 
is because when making the p -harmonic replacement of ft we gained new 
points of local injectivity for /i+i. These are all points in Ui. At the same 
time we did not loose local injectivity at the points where fi was already 
injective. 

• /iv+i : ^ Q is a homeomorphism, because it is a local homeomor- 
phism and the preimage of any point is a continuum in X. 


• For each z € U, we have \fN+i{z) — h{z) \ < 3e. Indeed, by triangle 
inequality, 


N 

\fN+i{z) - h{z) I < \fz+i{z) - fi{z) I 

i=l 


Let us take a quick look at each term |/i+i(z) — fi{z)\. If z £ Ui then 
both /j+i and fi lie in Qi ; therefore, \fi+i{z) - fi{z) \ ^ diam Qi < e. 
This term vanishes if z ^ Ui. But, for a given point 2 : there can be at 
most three pre-cells containing z. In other words the above sum consist of 
at most three nonzero terms, each of which does not exceed e . 

Step 5. Letting e small. We are now ready to proceed to the final construction 
of the mappings hj : X Y, 


( fN+i{z) , where N = N{e), e = l/j ii z £U 

[ h{z) if z £X\U 

Obviously, we have \hj{z) — h{z)\ ^ 3/j everywhere in X. Hence hj =| h 
uniformly in X. To complete the proof of Proposition 4.1, we need only verify that 
hj —>■ h strongly in ). The crucial ingredient is that <Sx[hj] ^ (Sx[h ], 

for all j = 1,2,... . In particular, hj converge to /i weakly in . Now 

the lower semicontinuity of the energy functional yields (ox[h] ^ liminf ^ 

liminf (S%\h] = S^fh] . This in turn implies, by uniform convexity arguments, that 
hj converge to/i strongly in 1^^’P(X, K^). □ 

Remark 4.4. Here is a careful look at the uniform convexity arguments. Consider 
the coordinate functions for hj =Uj+iVj and h = u+iv . In view of lower semicon¬ 
tinuity, we have f |Vm|p ^ liminf f |Vuj|P and f |Vn|^ ^ liminf f |Vuj|p . Adding 
these inequalities, we obtain f (|Vu|^ -I- |Vn|^’) ^ liminf f (|Vuj|P -I- |Vuj|P) = 
liminf ^[hj] ^ £’[h] = f (|Vm|p -I- |Vu|p) , which is possible only when f |Vu|^ = 
liminf and / |Vn|^ = liminf Now we see that Vuj and Vvj 

converge strongly in J§fP(X, M^), because the usual normed space .5f^’(X, R^) is 
uniformly convex, by Clarkson’s Inequality for vectors in the Euclidean space , 
see [12]. Actually, one could apply Theorem 2 in [13] to infer that the Banach 
space L^(X,R^) x L^(X,R^), equipped with the norm [/^d/j^ + , is 

uniformly convex as well. 


5. Completing the proof of Theorem 1.3 

We now return to the surfaces X G ^ and Y C and construct the mappings 
hj : X Y stated in Theorem 1.3. The arguments are similar to those used in 
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Steps 3,4 and 5. The main difference, however, is that we now choose and fix one 
particular finite cover of Y by cells. Having Proposition 4.1 in hands there will 
be no need to partition those cell into smaller cells. We adopt analogous notation. 
Thus, we let Y C be covered by Lipschitz cells Qi, Q 2 , ■■■, Qn , including both 
internal and boundary cells. This time N is fixed for the rest of our proof. In 
symbols, 

• Y= QiUQ2U...UQAr 

• 1 ^ E*=i Xq, (y) , for y G Y 

Let £ be any positive number. As before, we proceed by induction to define a 
chain Ai T 2 Tv F^+i of monotone mappings Fi : X Y. 

dcf 

The first map is = h. In the induction step we appeal to Proposition 4.1. 
Suppose we are given a monotone map Fi : X Y of Sobolev class #^^’P(X, Y) 
and its monotone extension Fi : ‘3/. Consider the cell Qi C Y and the 

corresponding pre-cell Ui G ^ for Fi. Then, by Proposition 4.1, there exists a 
monotone map Ti+i : X Y such that: 

(a) F,+i G “^(X, Y) n iri’P(X,]R3) 

(b) = A, : X\Wr-^ Y\Q,, 

(c) Fi+i : Ui Qi is a homeomorphism, 

(d) \F,+i - 2 ^ everywhere in X, 

(e) II DFi+i - DFi II^P(x) ^ 5 ^ 

In each induction step, passing from Fi to Ti+i we gain injectivity of Ti+i 
within Ui , and at the same time produce no branch points in the previous pre¬ 
cells. Thus Fm +1 '■ X Y is a local homeomorphism. Arguing as in Step 4, since 
Fiv+i : X Y is monotone, we see that Fm+i '■ X Y is a homeomorphism. 
Lastly, by triangle inequality, we obtain 

Hpf 

IIIPa-i-i -/i|||5?p(x) = II-Pa-i-i - T’i||<g>(x) + II-DFat+i - HT’ill^p(x) < 

N N 

I-Pi-I-I “ ^ill'^(x) + X/ “ DFi\^p(x) ^ 2 2 ^ ^ 

i=l i=l 

The proof of Theorem 1.3 is concluded by setting e = 1/j and hj = F^+i ■ 

6. Applications to Thin Plates and Films 


Let us demonstrate the utility of Theorem 1.3 by establishing the existence of 
the energy-minimal deformations of thin plates (planar domains) and films (sur¬ 
faces) for p -harmonic type energy. Recall the reference manifolds and the Jordan 
domains X C ^ and Y C . Here we assume that both X and Y are Lip- 
shitz domains. We examine homeomorphisms h : X —>■ Y in the Sobolev space 
C 1^^’^(X, Y), 2 ^ p < 00 , and their weak limits. Every homeo¬ 
morphism h G /^^’^(X, Y) extends up to the boundary as a continuous monotone 
map, still denoted by ft. : X Y. Moreover, we have a uniform bound of the 
modulus of continuity in terms of the Dirichlet and the (S’p- energy on a surface, 
see Subsection 6.1 below for the definition of £’p . 


( 6 . 1 ) 


|ft(a:i) - h{x 2 )\'^ 


CiX,Y)^2[h] 
log (e + ) 


CpiX,Y) [<S’p[h]f^^ 
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For h fixed, its logarithmic modulus of continuity was already known by Lebesgue 
[37]. However, it is the dependence on the energy of h that we are specifically 
concerned (to apply limiting arguments). The proof of estimate (6.1) runs along 
similar lines as for Lipschitz planar domains in [28]. There are, however, routine 
adjustments necessary to fit the arguments to 2 -dimensional surfaces. 

Now a -weakly converging sequence of homeomorphisms between X and Y 

actually converges uniformly. Based on Theorem 1.3, we have 

(6.2) :^(X, Y) = i^(X, Y) = Y) 

meaning that (respectively) the strong closure, the sequential weak closure, and the 
monotone maps in the Sobolev space /^^’^(X, Y) are the same thing. 

6.1. The isotropic p -harmonic integral on surfaces. An intrinsic example 
of the energy-minimal deformations of thin plates and films is furnished by the 
p -harmonic integral, 2 ^ p < oo. Suppose we are given a monotone Sobolev 
map h G ./#p(X, Y). To almost every point a; G X there corresponds the linear 
tangent map Dh{x) : Ta;(X) Ty(^), y = h{x) G Y, and its adjoint D*h{x) : 
Ty(^^) Ta;(X) with respect to the scalar products in Ta;(X) and Ty{‘3 ^). The 

Cauchy-Green stress tensor Gh = [D*h] o [Dh] : Ta;(X) —>■ Ta;(X) gives rise to 
the Hilbert-Schmidt norm of the tangent map, \Dh\ = [Trace G/j ]^/^ . Now the 
isotropic p -harmonic energy of h is defined by: 

(6.3) 'S’p[h] = J \Dh{x)\Pdx, 

where the area element da; is the one induced by the Riemannian metric in X. The 
term isotropic refers to the fact that the integrand is invariant under the rotations 
in Ta;(X) and Ty(Y). We call S’ 2 [h] the Dirichlet energy. The energy in (6.3) fits 
to the following more general scheme: 

(6.4) S'[h] = J Fi{x,h,Dh)dx , for mappings /i G ./#p(X, Y), 

where E(a;, y, L) is a given real-valued function defined for a; G .1?^ and y € '3^, and 
the linear maps L : Tx.S^ Ty3^ ■ We shall impose the following conditions on 
the energy integral in (6.4); they suffice for the application of the Direct Method 
in the Calculus of Variations: 

• coercivity; 

c/^ |Dh(a;)|Pda; ^ J^E{x,h, Dh) dx ^ |Dh(a;)|Pda;. 

• continuity in the strong topology of #'^’^’(X, Y) ; 

(S’ [h] = lim inf (S’ [hj] , whenever hj G ^p{X, Y) converge strongly to h. 

• lower semicontinuity, 

S [h] ^ lim inf (S [hj] , whenever hj G .^p(X, Y) converge weakly to h . 

These conditions hold, in particular, for the p -harmonic integral in (6.3). Next, 
choose and fix a homemorphism p G J^(X, Y) and a compact set F C 9X (empty 
in case of traction free problems). We consider the class of Sobolev homeomor¬ 
phisms, 

(6.5) J^p{X,Y,r;p) = {h £ JSfp{X,Y)-, hy = py (upon extension to i9X) } 
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By Theorem 1.3, its strong and sequential weak closures in are the 

same and coincide with the monotone mappings of the class: 

( 6 . 6 ) ^p{X,Y,r;^) = {h e ^p{X,Y); hir = <f\r } 


Straightforward application of the direct method in the Calculus of Variations yields 

Theorem 6.1. There always exists the energy-minimal map ho € ^p(X,Y,r; p) 
such that 


S [/lo] = min S \h] = inf \h\ (no Lavrentiev phenomenon) 

h£^p{x,Y,r-,(p) hGjiep{x,Y,r-,ip) 


Remark 6.2. The existence of energy-minimal mappings within ^p{X,Y), p ^ 
2, can also be obtained for many realistic variational integrals in NE, including 
neohookean energy; 

£[h] = ^ [\Dh\^ + for h € ^+{X,Y) 

Here ^2 stands for mappings in ^ 2 {X,Y) with positive Jacobian determinant. 
For ^-connected, 2 ^ £ < 00 , Lipschitz domains X and Y the class ^ 2 {X,Y), as 
oppose to J^(X, Y), is closed under weak convergence in 1^^’^(X, Y). Applying 
the direct method in the calculus of variations we see that there always exists an 
energy-minimal map ho G ^^(X,Y). Now, by Theorem 1.3 the minimizer ho can 
be approximated uniformly and strongly in by homeomorphisms. The 

key issue is whether 


(6.7) 


£[ho] = inf £[h] 
hGJif2(X,Y) 


It is interesting to notice that (6.7) holds if 

1 


det Dho 


G -^l“c(X). 


Indeed, this condition implies that ho has locally integrable distortion; that is, 
\Dho{x)\‘^ < K{x) detDho{x) , where K G .ifio(.(X). 

By [32] ho is open and discrete and, being monotone, is a homeomorphism of class 

jr2(X,Y). 


7. Afterward 


Remarkably, the existence of traction free minimal deformations requires ad¬ 
vanced topological arguments, as compared with the classical Dirichlet boundary 
value problems. Before further reflections let us take a brief look at one more ex¬ 
ample. 

An Example. Given a pair of circular annuli X = {z; r < jzj < i?} and 
Y = {w; 1 < |u;| < ^ {R R~^)} , r < 1 < R, the energy-minimal deforma¬ 
tion in ^ 2 {X.,Y), unique up to the rotations, is given by: 


h{z) = 



r < jzj ^ 1 

1 < l^l < i? 


( critical harmonic Nitsche map ) 
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see [3, 29]. Squeezing of {z; r<| 2 ;|<l}cX into the inner circle of Y manifests 
itself. Let us confront it with the solution of the classical Dirichlet problem under 
the same boundary values as h ; that is, 


(7.1) 


h{z) = 




The solution is a harmonic mapping 


for jzj = r 
for \z\ = R 


h{z) = A z + — 
z 


where A = 


R^-2r + l 
2{R^ - r2) 


and B = 


r{2R? — rR^ — r) 
2(i?2-r2) 


which, as expected, fails to be injective. The point is that, under this harmonic 
solution the annulus X folds along the circle jzj = ■ 

It is axiomatic in NE that folding should not occur under hyperelastic deformations. 
Uniform limits of hyperelastic deformations do not exhibit foldings as well. 

It is a common struggle in mathematical models of Nonlinear Elasticity [2, 5, 6, 7, 
8, 9, 11] to establish existence of the energy-minimal deformations which comply 
with the principle of no interpenetration of matter. The phenomenon of squeezing a 
part of a domain, as described by monotone mappings, is inevitable. It is therefore 
reasonable to adopt Monotone Sobolev Mappings as legitimate deformations in the 
mathematical description of 2D -elasticity (cellular mappings in higher dimensions). 
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